On the contribution of twist-3 multi-gluon 
correlation functions to single transverse-spin 
asymmetry in SIDIS 

Hiroo Beppu^, Yuji Koike^, Kazuhiro Tanaka^ and Shinsuke Yoshida^ 

^ Graduate School of Science and Technology, Niigata University, Ikarashi, Niigata 950-2181, 
Japan 

^ Department of Physics, Niigata University, Ikarashi, Niigata 950-2181, Japan 
^ Department of Physics, Juntendo University, Inzai, Chiba 270-1695, Japan 

E-mail: tanakakOsakura. juntendo . ac . jp 

Abstract. We study the single spin asymmetry (SSA) induced by purely gluonic correlation 
inside a nucleon, in particular, by the three-gluon correlation functions in the transversely 
polarized nucleon, . This contribution is embodied as a twist-3 mechanism in the coUinear 
factorization framework and controls the SSA to be observed in the _D-meson production 
with large transverse-momentum in semi- inclusive DIS (SIDIS), ep^ eDX. We define the 
relevant three-gluon correlation functions in the nucleon, and determine their complete set 
at the twsit-3 level taking into account symmetry constraints in QCD. We derive the single- 
spin-dependent cross section for the D-meson production in SIDIS, taking into account all 
the relevant contributions at the twist-3 level. The result is obtained in a manifestly gauge- 
invariant form as the factorization formula in terms of the three-gluon correlation functions and 
reveals the five independent structures with respect to the dependence on the azimuthal angle 
for the produced D meson. We also demonstrate the remarkable relation between the twist-3 
single-spin-dependent cross section and twist-2 cross sections for the D-meson production, as a 
manifestation of universal structure behind the SSA in a variety of hard processes. 



1. Introduction 

The single transverse-spin asymmetry (SSA) in the semi-inclusive DIS (SIDIS), e{£)+p{p, S±) ^ 
e{l') + h{P}i) + X, arises as T-odd effect in the cross section for the scattering of transversely 
polarized nucleon with momentum p and spin S^, off unpolarized lepton with momentum 
observing a hadron h with momentum in the final state. Here, q = (. — I' , Q"^ = —q^, and 
Q ^ Aqcd- The SSA can be observed also in pp collisions, the pion production p^p — )■ irX [Ij, 
and the Drell-Yan and direct-7 production, p^p — >■ 'j^*'^X [2]. Similarly as these examples, the 
SSA in the SIDIS requires, (i) nonzero transverse-momentum Pfij_ originating from transverse 
motion of quark or gluon; (ii) nucleon helicity flip in the cut diagrams for the cross section, 
corresponding to the transverse polarization S±; and (iii) interaction beyond Born level to 
produce the interfering phase between the LHS and the RHS of the cut in those diagrams. 
When Pfi± <^ Q, all (i)-(iii) may be generated nonperturbatively from the T-odd, transverse- 
momentum-dependent (TMD) parton distribution/fragmentation functions such as the Sivers 
function (see [3ii4j). By contrast, for large Ph± ^ Aqcd, (i) should come from perturbative 



mechanism as the recoil from the hard (unobserved) final-state partons, while nonperturbative 
effects can participate in the other two, (ii) and (iii), allowing us to obtain large SSA. This 
is realized with the twist-3 distribution/fragmentation functions in the collinear-factorization 
framework. In we derived the corresponding factorization formula, in the leading-order 
(LO) perturbative QCD, for the SSA associated with the twist-3 quark-gluon correlation inside 
the nucleon, and provided a practical procedure to calculate the relevant partonic hard parts 
manifesting their gauge invariance at the twist-3 level. A similar twist-3 mechanism for the SSA 
can also be generated by the purely gluonic correlations inside the nucleon. In particular, the 
SSA to be observed in the production of the D meson with large transverse-momentum Pfi± is 
controlled by the corresponding gluonic contributions, shown in figure [TJ through the photon- 
gluon fusion to create a pair of c and c quarks, one of which fragments into the D meson. In this 
report we discuss those gluonic effects, taking into account the nonzero masses, rUc and m/^, for 
the charm quark and the D meson j6]; our SSA formula for the high Ph± D-meson production, 
ep^ eDX, differs from the corresponding result in the literature [7], and the origin of the 
discrepancy will be clarified. 




(a) 



(b) 



Figure 1. Generic diagrams 
for hadronic tensor of ep^ — > 
eDX induced by gluonic ef- 
fects in the nucleon. Each one 
is decomposed into nucleon 
matrix element (lower blob), 
D-meson matrix element (up- 
per blob), and partonic hard 
scattering part by the virtual 
photon through the photon- 
gluon fusion (middle blob). 



2. Three-gluon correlation functions in the transversely polarized nucleon 

Figure [IJa) gives rise to the twist-2 contribution for the unpolarized SIDIS, ep — )• eDX, and 
the corresponding cross section is expressed by the factorization formula associated with the 
unpolarized gluon-density distribution for the nucleon, 

G{x) = -^-fJ ^e*^-(p|Ff"(0)Fr(An)b), (1) 

as the nucleon matrix element of the gauge-invariant lightcone bilocal operator using the 
gluon field-strength tensors, F°"(An) = F^l^{\n)np, F^l^ = - dl^ A'^ + gfabcA'^A^; here, 

p^ = {p^,0,0±) is the nucleon momentum regarding as lightlike up to the twist-3 accuracy, 
= (0,n~,O_L) is another lightlike vector satisfying p ■ n = 1, and we have suppressed the 
gauge- link operator which connects the field strength tensors so as to ensure the gauge invariance. 
Likewise, the "three-gluon distribution" functions, relevant to the SSA arising from figure [T][b), 
are defined through the gauge-invariant correlation function of the three field-strength tensors, 

■^f!IL(^1'^2) =-9^'J ^ / ^e^^-^e^'^^— ^)(p5|F,'^"(0)F7"(^n)Fr(An)|p5) 

= ^d''''0^f'^x^,X2) - ^r''^N"f'^xuX2), (2) 

for the nucleon with the transverse spin vector = (0, 0, Sj_), normalized as S"^ = — = —1. 
This may be regarded as an extension of the quark-gluon correlation functions discussed in, 



e.g., [5], with xi^2 = ^1,2 • n denoting the relevant momentum fractions of the gluons shown in 
figure mb). Here, d"''"^ and /"'"^ are the usual symmetric and anti-symmetric structure constants 
of the color SU(3) group, and the constraints from hermiticity, invariance under the parity and 
time-reversal transformations, and the permutation symmetry among the participating gluon 
fields allow the further decomposition in terms of the real, Lorentz-scalar functions, 0{xi,X2) 
and N(xi,X2), which are both symmetric under the interchange xi -H- X2, such that [6] 

= 2iMN [0(xi,X2)5"^e^f'^^ + 0{-X2,xi - X2)/^e"f'^^ + 0{x2 - xi, -xO^^^e'^f""^] , (3) 

and similarly for N"'^'^{xi,X2), at the twist-3 accuracy; here, e'^^""^ = e'^^'^pxn^Sy. 0(xi,X2) 
and N{xi^X2) are defined as dimensionless, accompanying the nucleon mass Mjy in ([3]), and are 
associated, respectively, with the C-odd and C-even combinations of the three gluon operators in 
([2]), satisfying 0(xi, X2) = 0{—xi,—X2) and A^(xi,X2) = —N{—xi,—X2). Indeed, 0{xi,X2) and 
N{xi,X2) constitute a complete set for expressing gluonic correlation effects inside the nucleon at 
the twist-3 level. We note that the other two functions, in addition to 0{xi,X2) and A^(xi,X2), 
were introduced in [8] to constitute a basis of three-gluon functions at the twist-3 level, but this 
basis proves to be redundant, i.e., the two functions 0{xi,X2) and N{xi,X2) form a genuine 
complete set [6]. The authors of [£] also pointed out that there are only two independent pure 
gluonic functions at twist-3, in the study of the evolution equations for the twist-3 distributions. 

3. Collinear expansion and gauge invariance for the three-gluon correlation effects 

We calculate the hadronic tensor W^^{p, q, Ph), represented by figured) We are interested in the 
contribution in which the c and c quarks are created through the photon-gluon fusion process and 
one of them fragments into a D (D) meson, so that the corresponding fragmentation function 
D{z), with z being the relevant momentum fraction, is factorized from W^i, as 

/dz 
^D{z)Wf,uip,q,Pc) , (4) 

where the summation over the c and c quark contributions is implicit, and Pc is the momentum 
of the c (or c) quark with mass mc, using the lightlike vector w of 0{1/Q) satisfying P^- w = 1, 
we have, p^ = Pj^ / z + rw'^, with r = {rn^z — m|/z)/2 to satisfy p^ = ml, and, at the leading 
twist-2 accuracy for the quark- fragmentation process, we set = p^ / {Ph'P)- We also note that, 
in the LO with respect to the QCD coupling constant for the partonic hard scattering parts, 
the contribution of figure [U^a) can not give rise to the SSA, due to the symmetry properties of 
the correlation functions of two gluon fields in the polarized nucleon (see Appendix A in [6|). 
So we shall focus on the analysis of figure mb). We work in Feynman gauge and apply the 
collinear expansion to the hard scattering part in figure [T][|b), keeping all the terms contributing 
up to the relevant accuracy of twist-3. We perform the corresponding calculation on the basis of 
generalization of our previous formulation [3] that was developed for the case of the SSA for the 
pion production induced by the quark-gluon correlation inside the nucleon. For simplicity, we 
shall henceforth write w^i, {p,q,Pc) of (HD as tt; {p,q,Pc), omitting the indices for virtual photon. 
The contribution from figure [TJb) to w{p,q,pc) can be written as 

w{p,q,Pc) = J^J ^,Sfj^,{h,k2,q,Pc)M^,l\h,k2), (5) 

where S'^^J')^{ki,k2,q,Pc) is the partonic hard scattering part (middle blob of figure mb)), and 

Mj/(fci,fc2) = g [ d^C [ d^7^e^'^ie^^'^-'^^^(j,S\Aimciv)A^ampS) (6) 



is the corresponding nucleon matrix element (lower blob), including one QCD coupling 
constant g. In ([5]), a real contribution relevant to the spin-dependent cross section occurs from 
an imaginary part of the color-projected hard part C'^™5'^|J,a(^1) k2,q,Pc) with C'^°' = if^'^"' and 
(jbca _ ^bca^ gince Cj^"'M'^^^{ki, ^2) are pure imaginary quantities [6]. This means that only the 
pole contribution produced by an internal propagator in the hard part can give rise to SSA. 




Figure 2. Feynman diagrams for the par- 
tonic hard part in figure [Hb), representing 
the photon-gluon fusion subprocesses that 
give rise to the "surviving" pole contribution 
for ep^ — 7> eDX in the LO with respect to the 
QCD coupling constant. The short bar on the 
internal c-quark line indicates that the pole 
part is to be taken from that propagator. In 
the text, momenta are assigned as shown in 
the upper-left diagram, where pc denotes the 
momentum of the c-quark fragmenting into 
the D-meson in the final state. The mirror 
diagrams also contribute. 



In the LO with respect to the QCD coupling constant, we find that the diagrams shown in 
figure [21 together with their mirror diagrams, give rise to the "surviving" pole contributions; 
here, the quark propagator with a short bar produces the corresponding pole contribution. The 
other pole contributions turn out to cancel among themselves after summing the contributions 
of all the leading-order diagrams for ([5]). With the assignment of the momenta ki and k2 of 
gluons as shown in figure [21 the condition for those poles is given by {pc — k2 + ki)^ — = 0. 
After we perform the collinear expansion and reach the collinear limit, ki xip (i = 1, 2), 
with xi,X2 and X2 — xi representing the longitudinal momentum fractions of the relevant three 
gluons, this condition reduces to xi = X2 and hence a pole of such type is referred to as the 
soft-gluon pole (SGP). In the following, we assume that S^^J,*)^ {ki, k2,q,Pc) in © represents the 
sum of the contributions of the diagrams in figure [2] and their mirror diagrams, in which the 
barred propagator is replaced by its pole contribution. We write the hard part S'^'l\ {ki, k2,q,Pc) 
as Sfj,ux (fci, /C2), for simplicity, suppressing the color indices a, 6, c and the momenta q and pc, 
and, correspondingly, we write M^^'^'^(fei, ^2) as M'^^^{ki, k2)- 

To perform the collinear expansion, as usual, in the "hadron frame" [5] with (7^ = {q^,q) = 
(0, 0,0, — Q), where both the virtual photon and the initial nucleon move along the z-axis, we 
decompose the relevant gluon momenta ki {i = 1,2) as /cf = {ki ■ n)p^ + {ki ■ p)n^ + kj_ = 
Xip^ + uj^j^ki, with u)^^j^ = g^j^ — p^^n^: Since p^ ~ g'i^Q in the hadron frame and thus 
the component along p^ gives the leading contribution with respect to the hard scale Q, 
we expand S^u\{ki, k2) around ki = Xip. Expressing also the gluon field A" as A" = 
(p^n^ + w"^) A'^ = p°'n-A + u}°'i^A'^, we note that, in the matrix element M^^^^{ki, ^2) of ([H), the 
components associated with the second term a;"^^'^ give rise to the contributions suppressed by 
~ \/Q or more, compared with the corresponding contribution due to the first term, p'^n ■ A 
(see, e.g., [5]). We organize the integrand of ([5]) according to the order counting based on 
those decompositions, keeping the terms necessary in the twist-3 accuracy. For this purpose, 
we need the Taylor expansion of the hard partonic subprocesses up to the terms with the 
third-order derivative, as S^u\{ki,k2) ^ S^ux{xi,X2) + uj°'^k1 dS^u\{ki,k2)/dkf\j^_^^^ H h 



{l/Q)u^,k^uj^,klu\k'i d^S^,x{ki,k2)/dk^dk^dk^^ _ , where 5^,A {xi,X2) = Sfj,ux{xip,X2p), 

because the terms with the third derivatives could produce the contributions that behave 
as the same order as the term, ^^^^(xi, X2)w"^a;'^^a;'^o.M'^'^'^(/ci, A;2), arising in the colhnear 
expansion of the integrand of ([5|): If we substitute Sai3-y{xi, X2)io°'i^uj^T-io'^f^M'^'^'^ {ki, k2) directly 
into the integrand of ([5|) and perform the integrals over ki and k2, the result produces the 
contribution, which is associated with the hard scattering between the three physical gluons 
from the nucleon and behaves as the same order as the formal convolution of Sap'y{xi,X2) with 
the twist-3 correlation functions in ([2]); this corresponds to a quantity of twist-3. Thus, our 
collinear expansion produces lots of terms, each of which is not gauge invariant, and it would 
look hopeless to reorganize those into a form of the convolution with only the gauge-invariant 
correlation functions 0"^'^(xi,X2) and N°'^'^{xi,X2) of ([2]) used. 

However, as was the case for the pion production associated with the twist-3 quark-gluon 
correlation functions [5] , great simplification occurs due to Ward identities satisfied by the 
corresponding partonic hard-scattering function, S^u\{ki, k2). To derive the Ward identities, 
we note that the contribution to S^px{ki, k2) from each diagram in figure [2] all contains the two 
delta functions, 6 {{k2 + q — Pc)"^ — ttiI) and 5 {{pc + ki — ^2)^ — m^c), representing the on-shell 
conditions associated, respectively, with the final-state cut on the unobserved c-quark line and 
with the unpinched pole contribution of the barred propagator. Also, in each diagram in figureEl 
the c-quark line fragmenting into the final-state D-meson is on-shell (see the discussion below 
dl])). Due to these on-shell conditions for the diagrams in figure [2] and the similar conditions for 
their mirror diagrams, S^y\{ki,k2) satisfies the Ward identities [U], 

fci5'^vA(^i,^2) = 0, k2Sf,yx{ki,k2) = {k2 - ki)^Sf,^x{h,k2) = 0, (7) 

which can be used to reorganize our collinear expansion. Indeed, combined with another relation. 



dS^,x{ki,k2)p' 



dkl^ 



dS^,x{kiM)p' 



ki=XiP ^-L 



(8) 



which can be derived by direct inspection of the diagrams in figure [2] and their mirror diagrams, 
we find [6] that all the gauge-noninvariant terms in our collinear expansion vanish or cancel 
among themselves, and the remaining terms can be expressed in a gauge-invariant form, as the 
factorization formula in terms of the three-gluon correlation functions of twist-3, defined in ([2]), 



w {p, q, Pc 



dxi f dx2 dS^l\{ki,k2,q,Pc)p 



xi J X2 dkn 



<a;^^u;'^,A1^5>i,X2), (9) 



up to the higher-order corrections beyond the present accuracy. Here, we have restored the color 
indices a, b, c as well as momentum variables q, pc, which were associated with the hard-scattering 
function in ([5]). Thus, the total twist-3 contribution to w{p, q,Pc), relevant to SSA, proves to be 
expressed solely in terms of the three-gluon functions 0{xi,X2) and A''(xi,X2) (see ([3])). 

When calculating the hard part in Q, dS'^'lf^{ki,k2,q,Pc)p^/dk2 , the derivative with 

respect to can hit the delta functions, 6 {{k2 + q — Pc)'^ — rn^) and 6 {{pc + ki — ^2)^ — w^), 
which are involved in S'^'l\{ki,k2,q,Pc)p'^ as mentioned above ([7|). Such derivatives of these 
delta functions can be reexpressed by the derivatives with respect to X2, and then be 
treated by integration by parts, giving rise to the derivative of the three-gluon functions 
0{xi,X2) and N{xi,X2)- After such manipulations, we have 6 {{pc + ki — k2 



\2 ^2 



m 



c) \ki=Xip 



{l/2pc ■ p)S{xi — X2), manifesting the SGP contribution, while another delta function becomes 
{x = Xi)j/x, z = Zf/z, qx = Ph±/zf with Xhj = Q'^/{2p ■ q) and Zf = p ■ Ph/p • q, as usual) 



6 {{k2 + 




k2=xp zQ"^ \x ) \z ) z^Q 



2 

rric 



--l)(--l)+-27% ' (10) 



implying x > xi,j] the similar contributions arise also from the corresponding mirror diagrams. 
Therefore, evaluating the SGP contributions in the hard part in ([9]), the resulting SSA proves to 
receive contributions associated with 0(x, x), 0{x, 0), N{x, x), N{x, 0), and also their derivatives 
with respect to x. In particular, the partonic hard part associated with 0{x,x) {dO{x, x)/dx) is 
different from the hard part associated with O(x,0) {dO{x,0)/dx), due to the difference in the 
tensor structures among the three terms in the RHS of Similarly, the partonic hard part 
for N{x,x) {dN{x,x)/dx) is different from that for N{x,0) {dN (x , 0) / dx) . 

Here, we make a brief comment on the calculation presented in [7] for the same phenomenon, 
i.e., for the twist-3 mechanism to the SSA in SIDIS, ep^ — )• eDX. The calculation by Kang and 
Qiu in [7j started from a factorization formula that was proposed as a straightforward extension 
of the corresponding factorization formula for the SSA in the pion production, associated with 
the twist-3 quark-gluon correlation functions, but Kang-Qiu's factorization formula did not 
manifest gauge invariance, nor permutation symmetry among the participating gluons (see 
equations (26) and (27) in [7j). Assuming gauge invariance in their factorization formula, Kang 
and Qiu claimed that w {p, q,Pc) can be eventually calculated with the following formula, 

^%9a/3>I^Jl(xi,X2), (11) 

in the notation used in this report, where g'^ = g^'^ — p^^n^ —p'^n^ = —S^S^ — ^t^P^s^upnS ^ ^^^^ 
we can make the replacement uj^^ — > up to the irrelevant corrections of twist-4 and higher. 

Combined with the property, S^gajs-M.'^^bc^^'^i^'^) ~ ^' "^l^ich follows from we see that 

the three-gluon correlation functions involved in (|lip are expressed by the two types of functions 
of X, Tq^\x, x), and their derivative, after evaluating the SGP at xi = X2 = x, where 

T^^\x,x) = J ^^^2Jll^e^^P^y^ ^gp^esynp{pS\C'r^^^ (12) 

with C^'^ = if^'^"', C^"" = d^'^"'; in contrast to our result mentioned above, the contributions 

Tq^\x,0), dTQ^\x , 0) / dx , associated with unequal arguments, do not arise. The functions 
are given by the contraction of ^ with the particular tensor g/sa^s-ynp- Using we have 

^T^+\x,x) = -4Mn [N{x,x) - N{x,0)] , ^Tif\x,x) = -AMn [0(x, x) + 0(x, 0)] . (13) 

Thus, the twist-3 SSA obtained in [7] implies the same partonic hard parts for 0{x, x) and 0(x, 0) 
(for dO(x, x)/dx and dO{x,0)/dx), and similarly for iV(x, x) and— A^(x,0) (for dN{x,x)/dx and 
—dN(x,0)/dx); clearly, such result contradicts with the above-mentioned result based on our 
complete formula ([9]). It is straightforward to see that, if the tensor structure of the three- 
gluon correlation function A^^^£j'[c(xi, X2) of were assumed to be given by only one structure, 
g<^l3^ipnS ^ our formula ^ would reduce to the formula (jlip . up to the corrections of twist-4 and 
higher. However, such assumption contradicts with the permutation symmetry required by the 
Bose statistics of the gluon, as discussed in section 2 and represented in 

4. Result for the tvi^ist-S single-spin-dependent cross section for ep^ eDX 

We calculate the hadronic tensor @ using our factorization formula ([9]), and evaluate the 
contraction L^^Wfj,^ with the leptonic tensor for the unpolarized electron, L^j^y = 2{i^£'^^+£yi'^) — 
Q^9fj,u- The result is decomposed according to the dependence on the azimuthal angles (/) and $5 
for the initial-lepton's 3-momentum i and the nucleon's spin vector S, respectively, measured 



dxi 

Xl 



dx2 dSfilih,h,q,Pc)p^gT 



X2 



from the axis along the transverse-momentum Pfi± in the hadron frame, and, through standard 
manipulations, leads to the spin-dependent, differential cross section for ep^ — )• eDX [6]: 
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dO{x,x) 20{x,x) 



dx X 
dN{x,x) 2N{x,x 



dO{x,0) 2O(x,0) 



dx 



a4 



dx X 
dN{x,0) _ 2N{x,0) 
dx X 



2 , 0(X,X)^^3 _^ O(x,0) 



A<7^ + 



. ^2 N(x,x) . „o 

X 



Aa^ ,(14) 



where Sep = we use the shorthand notation, [da;] = dxhjdQ'^dzfdq^dcj), for the differential 

elements, and the summation X]fc implies that the subscript k runs over 1,2,3,4,8,9, with 



Ai = 1 + cosh^ if}^ A2 = —2, A3 = — cos <f> sinh 2^/;, A4 = cos 2(/» sinh^ ip, 



As = — sin (f) sinh 2^/;, ^Ig = sin 2(f) sinh^ '(/', 



(cosh Tp 



2xi)jS, 



1 



(15) 



and Sk defined as Sk = sin<I>5 for k = 1, 2, 3, 4 and Sk = cos $5 for k = 8, 9. The quark-flavor 
index a can, in principle, be c and c, with 6c = 1 and (5c = —1, so that the cross section for the 
Z)-meson production ep^ — )> eZ)X can be obtained by a simple replacement of the fragmentation 
function to that for the D meson, Da{z) — > Da{z). The delta function of ()10p appears, Oem is 
the fine-structure constant, = g^/(47r) is the strong coupling constant, and Cc = 2/3 is the 
electric charge of the c-quark. Partonic hard parts Aa^ depend on rric as well as other partonic 
variables; for the explicit formulae of Acj^, we refer the readers to [B]. We find that Acj^ 7^ Aa^, 
Acj^ 7^ A(j^, while the hard parts arising in the second and third lines in (|14p are common. 

The single-spin-dependent cross section p4|) can be decomposed into the five structure 
functions, based on the different dependences on the azimuthal angles $5 and (p through the 
above-mentioned explicit forms of Ak and Sk- The five independent azimuthal structures of this 
type have been observed also in the twist-3 single-spin-dependent cross section for ep^ — )• evrX, 
generated from the quark-gluon correlation functions, as presented in [10^ lllj. Introducing the 
azimuthal angles (ph and (ps of the hadron plane and the nucleon's spin vector 5, respectively, 
as measured from the lepton plane, they are connected to the above and (p as = (ph — (ps^ 
(p = (ph, and one may express ([T3|) as the superposition of five sine modulations, 

d^A(T 



[du] 



fi sm{(ph - (ps) + /2 sm{2(ph -(ps) + f3 sin (ps + sin(3(/),, - (ps) + h sin(</>/i + (ps), (16) 



with the corresponding structure functions /i, /2, • • • , /s- A similar form of five independent 
azimuthal dependences in the SSA was obtained in the TMD approach [3j for small- regions. 

It is worth comparing p4|) with the unpolarized cross section for SIDIS, ep — >■ eDX, which can 
be obtained straightforwardly by applying the above-developed formalism to extract the twist-2 
contribution from figure [TJ^a), generated from the unpolarized gluon-density distribution ([T]), as 



^S^unpol 

\duj] 




1-^1 



mt. 



\J2 Da{z)Gix)al! 



(17) 

in the same notation as in (|14p . with = = 0. Comparing the explicit form [6] of the LO 
partonic hard cross sections with the partonic hard parts Ao"^ in p4p . associated with the 
derivatives, dO{x,x)/dx and dN{x, x)/dx, of the three-gluon functions, one finds the following 
relations between the partonic hard parts at the twist-3 level and those at the twist-2 level: 



^4 



:^k ■ 



Q2(l_^) 



(18) 



The other partonic cross sections at the twist-3 level, A(t| (i = 2, 3, 4) in ([Tl]) . are also related to 
the partonic hard scattering parts at the twist-2 level, although, unlike (jlSp . the corresponding 
relations cannot be manifested by direct comparison between their explicit formulae. We find 
that all these remarkable relations are consequences of the "master formula" , 

which allows us to calculate ([9]) using the partonic hard part T-L'j^lxp, q,Pc) for the 2 — > 2 Born 
subprocess, which coincides with the LO contribution to the hard part in figure [U^a), up to the 
color structure. The proof of ()19p will be presented elsewhere. The similar master formula was 
derived \W\ [12] for various processes associated with twist-3 quark-gluon correlation functions. 

5. Summary 

We have investigated the SSA for the D-meson production in SIDIS, generated from the 
twist-3 three- gluon correlation functions for the nucleon. We first showed that there are only 
two independent three-gluon correlation functions of twist-3, 0{xi,X2) and N{xi,X2), which 
correspond to two possible ways to construct color-singlet combination composed of three 
active gluons. Then, we have formulated the method for calculating the twist-3 single-spin- 
dependent cross section generated from the three-gluon correlations. Our formulation is based 
on a systematic analysis of the relevant diagrams in the Feynman gauge and gives all the 
contribution to the cross section at the twist-3 level in the LO in perturbative QCD, guaranteeing 
the gauge invariance of the result. As a result, the SSA occurs as the pole contribution of an 
internal propagator in the partonic hard-scattering subprocess, and the corresponding SGP 
contribution leads to the SSA expressed in terms of the four types of gluonic functions of the 
relevant momentum fraction x: 0{x,x), O{x,0), N{x,x) and N{x,0). We find that all these 
four-types of functions and their derivatives with respect to x contribute to the final form of the 
single-spin-dependent cross section, generating five independent structures on relevant azimuthal 
angles. We also mentioned the master formula that manifests universal structure behind the 
SGP contributions, allowing us to relate the 3^2 subprocess relevant for the twist-3 level to the 
2 — > 2 subprocess for the twist-2 level. Our formalism developed for the twist-3 mechanism of the 
SSA arising from multi-gluon correlations, and also the above features revealed in ep^ — )■ eDX, 
apply to the SSA in the other processes, like An in p^p hX {h = vr, K, D, etc.) [T3|. 
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